Zero bias anomaly out of equilibrium 
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The non-equilibrium zero bias anomaly (ZBA) in the tunneling density of states of a diffusive 
metallic film is studied. An effective action describing virtual fluctuations out-of-equilibrium is 
derived. The singular behavior of the equilibrium ZBA is smoothed out by real processes of inelastic 
scattering. 
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The suppression of tunneling current at low bias due 
to electron-electron interaction is known as the zero bias 
anomaly (ZBA). The theory of ZBA for disordered metals 
at thermal equilibrium has been developed, on a pertur- 
bative level, by Altshuler and Aronov [l|, 0]- The non- 
perturbative generalization of this theory was achieved 
by Finkelstein Q . Measurements of the tunneling density 
of states (DOS) in biased quasi-one-dimensional wires |4[ 
call for an extension of the theory to non- equilibrium 
setups. In this work we study the ZBA for disordered 
metallic films out of equilibrium, in both the perturba- 
tive and the non-perturbative (in interaction) limits. 

Besides the experimental motivation, the problem of 
ZBA in a non-equilibrium system is of fundamental the- 
oretical interest. At equilibrium, the distribution of elec- 
trons in phase space has a single edge at the Fermi sur- 
face. The Coulomb interaction between the tunneling 
electron and the electrons in the Fermi sea excites vir- 
tual particle-hole pairs around the Fermi edge, leading 
to the suppression of the tunneling DOS, similarly to the 
Debye- Waller factor. The suppression gets stronger when 
the electron energy approaches the Fermi energy. Out of 
equilibrium, the distribution of particles may have sev- 
eral sharp edges rather than a single one at the Fermi 
surface, which poses important questions addressed in 
this work: How will the excitation of electron-hole pairs 
in this situation affect the tunneling DOS? Will there 
be an interpaly between the two edges? We show that 
the two edges are not independent: one edge affects the 
ZBA near the other via real interaction-induced scatter- 
ing processes governing the dephasing of electrons in the 
non-equlibrium regime. From this point of view the prob- 
lem we are considering is a representive of a class of phe- 
nomena that involve renormalization away from thermal 
equilibrium, such as the Fermi edge singularity [B[ and 
the Kondo effect @. 

What makes the ZBA particularly interesting is 
its deep connection to various conceptually impor- 
tant phenomenological ideas. At equilibrium, the non- 
perturbative results Q have been reproduced by quan- 
tum hydrodynamical methods Q, and, within the frame- 
work of the theory of dissipation Q, by methods that 
rely on the fluctuation-dissipation theorem. Our work 



circumvents this restriction. Starting with the Keldysh 
non-linear er-model, we derive an effective action that ac- 
counts for virtual fluctuations in disordered metals away 
from equilibrium. This action is complementary to the 
one for kinetics of real fluctuations (such as noise) devel- 
oped earlier HUME! . Il4| . Further, we discuss a connec- 
tion between our theory and phenomenological methods 
0> B B EH ■ As a central application of our theory, we an- 
alyze the ZBA problem and calculate the tunneling DOS 
for a two-dimensional (2D) diffusive metallic film subject 
to external bias. 

Consider an electron which tunnels from a tip into a 
metal, subject to an external bias. The Coulomb inter- 
action U(q) (U(q) — 2ne 2 /q, d = 2) causes electrons in 
the metal to readjust their position, such that they try 
to screen the added charge. Their motion is impeded by 
static disorder. For a sufficiently good metal, character- 
ized by the dimcnsionlcss conductance g = epT 1, (ep 
is the Fermi energy and r - the elastic mean free time) 
the electron motion is diffusive. Low energy excitations 
are accounted for by the non-linear Keldysh cr-model [l5[ 

iS[Q, 0] = fR^tr 1 + voM} 



-^Tr{D(VQ) 2 



4d t Q} - iTtv Tr{<p al a Q}, (1) 



where Q(r,t,t') is a 2 x 2 matrix field (in the Keldysh 
space) describing slow electronic modes, <j>(r, t) is a 2- 
vector in the Keldysh space representing the Coulomb 
interaction, and Tr denotes the trace both over Keldysh 
indices (Tr^ ), spatial, and temporal coordinates. The Q 
field is subject to a non-linear constraint, 



*iQ(r, i, h)Q(r, h,t') = 5{t - t>o, 



(2) 



7i = (To is a 2 x 2 unit matrix, 72 = o~\ the first Pauli ma- 
trix, is the DOS at the Fermi energy in the absence of 
interaction, and D is the diffusion coefficient. The theory 
can be rewritten in terms of the field Q only. Integrating 
the cf> field out, we derive the following action: 



= -^(Tr{I?(VQ) 2 -4d t Q}) 



iS[Q] = - 

2 Tr^{crig} [U 



L Tt k {Q}. (3) 



2 



Physical observables can be found by differentiating the 
generating functional 



Z[(pi,f2] = exp {iS[Q}) exp 



l + u U 

2<^i^2 - TTViTrif {Q} - n(p 2 Tr K {aiQ} 



(4) 



over the source fields ip± and ip2- 

So far all low energy excitation in the problem have 
been kept indiscriminately. The price one pays for this is 
an extreme technical complexity of the theory. In some 
cases, this description is excessive and the theory can be 
simplified by singling out a particular subspace of the 
Q matrices. The resulting theory is less general, but 
more suitable for tackling a specific class of problems. 
The noise statistics in disordered systems is a remark- 
able example of such a simplification. An extension of 
the Boltzman-Langevin approach (l6| to high order cor- 
relators is achieved on the subspace of matrices that are 
diagonal in the Wigner representation, 



Q(r,e,t) = e~ 



= '^f(r,e,t) A 2-4/(r,e,t)\ 



(r,M) 



(5) 



The theory llj, |l3|, [l4| resulting from the a- model ([3]) , 
reduced to the subspace ([5]), accounts for real fluctuations 
in agreement with the cascade idea of Nagaev 
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Below we use a similar ideology, "projecting" the er- 
model on the subspace appropriate for the ZBA prob- 
lem. According to works on ZBA at equilibrium, these 
are virtual fluctuations of gauge fields that dominate the 
suppression of the tunneling DOS BE El- This sug- 
gests that gauge-type fluctuations constituting local-in- 
time rotations of the saddle point A t _e, 



1 2-4n(r,t-i' 



A-t-i' = 







(6) 



where n(r, t — t') is the Fourier transform of the local dis- 
tribution function n(r, e), are to be retained. Plugging 
© into the action ([3]) and expanding up to quadratic 
order in / and /, we derive an effective theory of fluctu- 
ations in an interacting diffusive conductor. 



iS[fJ} = 



I'd 



(duj)(dq) 



f(-uj,-q)cj -Dq 



1 + voU(q) 



-Dq 2 T(uj)f(-oJ,-q)f(cu,q) 



(7) 



Here the effective temperature T(oj) is defined as 

T(w) = de n(e) [2 - n{e - uS) - n(e + w)] , (8) 



and n(e) is a distribution function determined by the 
Boltzmann equation. 

To demonstrate the consistency of our approach, we 
first show that the effective action (0 reproduces cor- 
rectly known density correlation functions. The sym- 
metrized density correlation function is given by 



1 1 d 2 Z 
o([pM),p(0,0)] = — — 77T - 7 rrrn \ 

2 J+ 4d0 2 (w,g)c'02(-w, -q) 



Similarly, the response function can be expressed as 

i(\p(T,t),fi(0,0)]-6{t)) u , q = ^-.jr— r 

2i 0<pi{<jj,q)0(p2{-u, —q) 



(9) 



2 V l + u U 



1 + VqU 



(10) 



Calculating the correlation functions entering Eqs. 
(fTU)) using the action ([7]), 



if f)uj,q 



ADq 2 T{uo) 



Dq 2 



\fu,qf—u), — q) — (ff)u,q 



l + u V(q) 
2 



to 



(Dq 2 



1+voU 



we reproduce the known results for for the spectral func- 
tion of density fluctuations 

1 r i voDq 2 T(uj) 

2 ( [/?(r ' t] ' ° )] = \D?[l + V0 U{ q )\+iu\* ' 

and the density-density response function, 



i(\p(r,t),p(0,0)]-m)u, q 



u Dq 2 



Dq 2 [l + v Q U{q)]-iuJ 



It is worth noting the analogy between our effective ac- 
tion and phenomenological theories that describe ZBA at 
equilibrium within an effective environment model 0, [Io| . 
The latter explain the ZBA as the influence of virtual 
fluctuations of an electromagnetic field at equilibrium [j| 
on the electron tunneling. In view of the fluctuation- 
dissipation theorem, fluctuations of the electromagnetic 
field are determined by the complex impedance of the 
system. In the zero dimensional case, modes of the elec- 
tromagnetic field can be considered as independent quan- 
tum harmonic oscillators. Being suddenly shaken by the 
incoming electron, they move away from the equilibrium 
position, reducing the overlap with their original config- 
uration hence suppressing electron tunneling amplitude. 
The action (TZJ) describes similar processes in a 2D diffu- 
sive system and without the assumption of thermal equi- 
librium, thus keeping the information about both real 
and virtual processes in a non-equilibrium state. 
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Now we are ready to apply the theory to the problem 
of ZBA out of equilibrium. The tunneling DOS 



"(e) = 2^ J(d P )(G R (e,p)-G A (e,p)) (12) 
can be rewritten in terms of Q matrices as 

v(z)=2j dt e~ ic (* _t ^(Trff Q(r, t', t)a z )Q . (13) 
Plugging Eq. (J6]) into Eq . (fT3)) . we find 
^ = l + 2i J dtn(t) e^-i'/'W sin fi Iff (t) \ , (14) 

where we used the notations 

/,/(*) = / {duj)(dq) (1 - COS^t) «//» W|9 , 



I ff (t) = / (du)(dq) sin tot ((ff)) 



(15) 



There is a subtle point related to the definition of the cor- 
relator ((//)) and ((//)) in JTSJ). The ZBA should vanish 
in the absence of interaction ({7 = 0). This physically ob- 
vious statement is valid to all orders of the diagrammatic 
calculations (all potentially singular contributions vanish 
after the time integration over the Keldysh contour) and 
is satisfied by the non-linear a model. However, it ceases 
to be an exact feature of the theory when the full space 
of Q matrices is reduced to the subspace Eq. (J6j) . This 
problem is easily cured. Since we are only interested in 
summing up the leading interaction-induced In 2 e contri- 
butions to the tunneling DOS (to all orders), these are 
the interaction-induced parts of (//) and (//) that we 
have to keep in Eq. |(TSJ) 0: 



(16) 



{(ff))u,q = (ff}u,q - (ff)u),q,U=Q 

-2iU 



Analogously, 



((//)) 



(Dq 2 - iu)[Dq 2 {\ + u U) - iw] 



ADq 2 T{uj){2 + is U)U 
\Dq 2 - iuj\ 2 \Dq 2 {l + v U) - iu\ 2 



(17) 



Equations (|14 p -(|17 p constitute our general result for the 
non-equilibrium ZBA. At equilibrium they reproduce the 
known results for the ZBA in the perturbative [l], Q and 
non-perturbative 0] regimes. 

To illustrate the effect of non-equilibrium conditions on 
ZBA we consider a diffusive film of length L connected to 
two zero-temperature reservoirs with voltage difference 
V. We assume that the energy relaxation in the film 
can be neglected, i.e. that L is shorter than the energy 
relaxation length, a condition met at not too high bias. 
In this case the solution of the Boltzmann equation is a 
double-step function 

n(e, x) = an (e - eV/2) + (1 - a)n (e + eV/2) , (18) 



where no{e) is the Fermi distribution function, a = x/L, 
and x is the distance from the reservoir. This distribution 
corresponds to the effective temperature 

T(w) = [a 2 + (1 - a) 2 ]T cq (iu) 

+a(l - a)[T cq (uj + eV) + T cq (uj - eV)} , (19) 



where T eq (w) is the equilibrium value 

T cq (aj) = lo coth(w/T) — > \u> 



(20) 



Substituting Eqs. (TH, (HH into Eqs.([l5},( EJ), we find 
the tunneling DOS out of equilibrium: 



V{e) ( 1 log-logpVrt- 

T 



- aex P \ a 2 n 



1 



+ (l- a )cxp^-^-^log^log(/J 2 K 4 rt + )J . (21) 

This result is valid in both the perturbative (when the 
exponentials can be expanded up to the first non-trivial 
term) and non-perturbative regimes. The DOS consists 
of two terms corresponding to electrons coming from the 
left and right reservoir. The energy scales governing the 
argument of the logarithm in these two terms are 



t^} = max 



<z(l-a)eVlog( 



eV 



,(22) 



respectively (k — 2ire 2 v). 

The evolution of the ZBA in the tunneling DOS with 
decreasing conductance (i.e. from the perturbative to 
the non-perturbative regime) is illustrated in Fig.[T] The 
DOS has a two-dip shape with minima reached at the 
energies where the distribution function exhibits discon- 
tinuous jumps (i.e. at the positions of the Fermi edge 
in the left and right leads, e = ±^-)- Away from the 
minima, the DOS is controlled by the energy measured 
from the corresponding edge. As this energy decreases 
(we get closer to one of the Fermi edges), the singularity 
near this edge gets affected by the presence of the other 
edge. The broadening of the ZBA singularity takes place 
on a new energy scale, 



, a(l^ /^V 

V ' A-kvqD S I eV 



(23) 



The notation introduced in Eq. (|23|) stresses the analogy 
with the equilibrium dephasing rate t^ 1 {T) governing 
the infrared cutoff of the interference phenomena such 
as the weak localization J2j. The emergence of t"^ 1 (V) 
shows that inelastic scattering processes (responsible for 
dephasing) lead to smearing of the ZBA singularity. 

Qualitatively the results can be explained as follows. 
Virtual fluctuations of the gauge fields act on a quasipar- 
ticle and suppress the tunneling DOS. Their influence is 
limited by the quasiparticle life time. An external bias 
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FIG. 1: Tunneling DOS vs. energy for the double-step dis- 
tribution function (|18[) with a — 0.5. The dimensionless con- 
ductance is, from top to bottom, g = 100, 10 (a) and 1 (b). 

enhances fluctuations of the electromagnetic field, which 
in turn dephase the tunneling quasiparticle. This results 
in a competition between the virtual processes establish- 
ing the ZBA and the real fluctuations (noise) cutting it 
off. 

Let us stress the crucial difference with the equilib- 
rium situation. At equilibrium, the ZBA singularity 
in 2D is cut off by temperature (which enters via the 
distribution function); specifically, the dcphasing rate 
r ( ^ 1 (T) does not affect the ZBA in any essential way, 
since t ( ^ 1 (T) <C T. It is thus a distinct feature of the 
strongly non-equilibrium regime that the ZBA is smeared 
by the inelastic scattering (dcphasing) rate. 

In conclusion, we have developed an effective theory 
for virtual fluctuations in disordered metals out of equi- 
librium. Using this theory, we studied the ZBA in a 2D 
metallic film biased by external voltage. We have found 
that out of equilibrium the tunneling DOS has a double- 
dip structure with minima reached at the " edges" of the 
particle distribution. The ZBA near any of the "edges" is 
influenced by the other one. The suppression of DOS is 
smoothed out by the real processes of inelastic scattering 
(dephasing) with characteristic energy scale T^iV). Fur- 
ther applications and extensions of our theory include, in 
particular, the ZBA in quasi-one-dimensional and strictly 



one-dimensional (Luttinger liquid) wires [20| out of equi- 
librium; these results will be reported elsewhere. 
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